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Abstract. We have undertaken a study of diluted magnetic semiconductors 
Gai-xMn^N and Gai^xCr^N with x = 0.0625, 0.125, using the all electron linearized 
augmented plane wave method (LAPW) for different configurations of Mn as well as Cr. 
We study four possible configurations of the impurity in the wurtzite GaN structure 
to predict energetically most favorable structure within the 32 atom supercell and 
conclude that the near-neighbor configuration has the lowest energy. We have also 
analyzed the ferro-magnetic as well as anti-ferromagnetic configurations of the impurity 
atoms. The density of states as well as bandstructure indicate half metallic state for 
all the systems. Tc has also been estimated for the above systems. 
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1. Introduction 



Gallium nitride is one of the most promising materials among the diluted magnetic 
semiconductor (DMS) material for application in spintronics. By doping transition 
metal (TM) atoms, Mn or Cr, local magnetic moment are introduced in semiconductor 
which mediate ferromagnetically. (Ga,Cr)N based DMS was predicted to show high 
Tc PP for high enough concentration of Cr and further Hashimoto et al. j2| observed that 
(Ga,Cr)N based DMS grown by ECR molecular beam epitaxy showed Tc above 400K. 
Cr "'"-implanted GaN, studied by photoluminescence and superconducting quantum 
interference device (SQUID) reveal that the implanted Cr"*" incorporates substitutionally 
at Ga site and the ferromagnetic order is retained upto 300K 1,3^. Takeuch et al. jl]. 
have reported a systematic study of changes in the occupied and unoccupied N-partial 
density of states (DOS) and confirm the wurtzite N 2—p DOS and substitutional doping 
of Cr into Ga sites using SXES and XAS. Recently, ferromagnetism above 900K was 
reported in Cr-GaN thin films Theoretically it was predicted that the ferromagnetic 
(FM) interaction in (Ga,Mn)N may be retained upto room temperature jHj. The initial 
reports of high Tc in (Ga,Mn)N were followed by controversial results where the reported 
Tc varied between 20K - 940K [3 IHl El QUI- Zajac and coworkers observed Mn ions in 
Gai_3;Mn2,N [x < 0.1) crystals coupled anti-ferromagnetically (AFM) JT]. Electronic 
structure and magnetic properties of zinc blende Gai_2^Mna;N for several values of x 
with varied spatial distribution of dopant atoms to understand the magnetic interaction 
for explanation of FM-AFM competition is discussed by Uspenskii et al. ^2] where the 
calculations were done using the tight binding LMTO method in the local spin density 
approximation. Sanyal and Mirbt ^H] have studied Mn doped GaAs and GaN DMS 
using the ab-initio plane wave code (VASP) within density functional theory (DFT). 
They have determined the interatomic exchange interactions by substituting Mn in 
various positions in the unit cell and have attributed the origin of ferromagnetism in 
(Ga,Mn)N to double-exchange mechanism involving the hopping of Mn— ci electrons. 
Raebiger et al. ^1] used the full potential linearized augmented plane wave (FP- 
LAPW) method to investigate the interplay between clustering and exchange coupling 
in magnetic semiconductor Gai-^jMua-As. They have studied all possible arrangements 
of the two Mn atoms on Ga sublattice for x ~ 6% and found that clustering of Mn 
atoms at near neighbour Ga sites is energetically preferred. Our analysis of the wurtzite 
GaN doped with Mn or Cr is motivated by the latter study. 

2. Method and Computational Details 

We have employed the spin-polarized Linearized Augmented Plane Wave Method (FP- 
LAPW) as implemented in the WIEN2K package J3] with the Generalized Gradient 
Approximation (GGA) for the exchange-correlation potential proposed by Perdew, 
Burke and Ernzerof (PBE96) ^H]- This is state-of-the-art electronic structure method, 
which does not use any shape approximation for the potential, to solve the Kohn-Sham 



type of equations self-consistently. 

GaN normally occurs in the wurtzite structure with lattice constants a = 3.19A and 
c = 5.19A, giving c/a ratio of 1.62. Each Ga is tetrahedrally bonded to N atoms at 
an average distance of 1.95 A and each N in turn is surrounded by four Ga neighbors. 
The calculations for DMS were performed within a 32 atom supercell, constructed from 
2x2x2 standard unit cell of wurtzite structure wherein the dopant is substituted at 
various cation sites, since it has been shown that the formation energy for interstitial 
Mn doping is higher than substitutional doping [T3]. The supercell approach is used to 
restrict the dopant concentration to a small value, which is of interest for studying 
magnetic properties of the system, without altering the original underlying lattice 
structure. Our interest was in observing the changes in the electronic structure of the 
DMS with respect to the possible different geometries of the dopants within the host 
semiconductor. Self-consistent electronic structure calculations were performed using 
the APW + local orbitals (lo) basis set for the valence and semi-core electrons with 
Rmt * Kmax = 7, Imax = 10 and total energy convergence of O.OOOOli??/. The Muffin-Tin 
radii for Ga, Mn and Cr were kept at 2.0A and that for N at 1.6A. Spin-polarized 
calculations were carried out to observe the effect of spin-splitting and to calculate the 
on-site magnetic moment at TM site. 

We have studied wurtzite GaN doped with one TM atom impurity, which is ~ 6.25% 
doping and two identical TM atoms in the 32 atom unit cell amounting to ~ 12.5% 
doping. To simulate different surroundings for the transition metal (TM) atoms we 
have spanned certain geometries within the 32 atom unit cell wherein the distance 
between the dopants is varied. In the case of the single impurity substitution, the 
nearest distance between two TM atoms is 6.38A in plane and 10.38A along the c-axes. 
We have studied four different geometries of two TM atom substitutions at 3.19A, near 
neighbor (nn), 4.5A , 5.19A and 6.08A separations. When two near neighbor (nn) Ga 
atoms are substituted by TM atoms the in-plane TM-TM atoms distance is 3.19A and 
along the c-axes it is 10.38A. For the second case the out-of-plane distance is 4.5A and 
the in-plane distance between the dopants is 6.3A . In the third case, two Ga atoms lying 
one above the other, along the z-axes, separated by a distance of 5.19A are substituted 
by TM atoms and the in-plane TM atoms are at 6.38A . The last case is such that 
the in-plane separation (6.38A) and out of plane separation (6.09A) between the TM 
atoms is comparable. For estimating the magnetically favorable system, the spins of the 
dopants are aligned along the same direction, corresponding to the FM configuration, 
and aligned in the opposite directions corresponding to AFM configuration. The self 
consistency was achieved on a mesh of 5 x 5 x 2 k-points. 

Structural relaxation for the TM site and the nn N sites was carried out to observe 
changes in the bond lengths between TM and the first shell of N atoms. Very small 
change (~ 3 — 4%) was observed in the bond lengths and no significant changes were 
seen in the band structure in agreement with the earlier reported results J7] . Thus the 
calculations reported here are for systems without allowing any relaxation. 



3. Results 



In the wurtzite GaN semiconductor, each Ga (N) is tetrahedrally bonded to 4 N (Ga) 
atoms. Pure GaN is a direct band gap semiconductor with top of the valence band 
consisting of N p-states and the bottom of the conduction band having Ga sp-character. 
The Ga d levels are deep and do not take part in the bonding. Thus they are treated 
as core states. The band gap of GaN, which is underestimated by density functional 
theory within the approximation used for the exchange correlation energy functional, 
is ~ 2eV. The experimentally determined band gap of undoped GaN is ~ 3.4eV. Das 
et al. ^H! have shown that, for Mn atoms to couple ferromagnetically, they need to be 
kept apart by more than the critical distance of 2.5A. Similar calculations on clusters 
of (GaN)Cr indicate that the critical Cr-Cr distance is 2.7A ^H]- In ah our calculations 
the distance between the dopants was greater than the corresponding critical distances. 



3.1. Mn doped systems 

Localized magnetic moments are introduced within the GaN system by substituting the 
cations with TM impurity atom(s). Mn atom with 4s^ and 3d^ electrons in the valence 
region replaces Ga atom with valency 4s^4p^. On substitution Mn atoms contributing 
five d levels per atom are thus expected to contribute to the observed magnetic moment. 
Since three of the valence electrons from Mn go into compensating the three electron 
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Figure 1. Mn-projected majority spin (a) dt^ and (b) bands in Gai^MnNiQ. 



states of substituted Ga, one hole per Mn is introduced into the system. Figure [T] shows 
the Mn-projected dt^ and dg majority spin electronic structure for Gan^MnNiQ. The Mn- 
d states lie at the top of the valence band and cross the Ci;' in some places. These are 
split into de and dt2 states, the dt2 level is two thirds filled and the de is almost occupied. 



The minority spin levels are empty and lie above the ep indicating 100% spin-polarized 
states. The Mn induced states lie in the gap region of GaN. The top of the valence band 
in GaN is composed of the N-p levels and the unique properties, particularly the half 
metallic state of DMS, thus arise from the TM d and host p interactions that couple the 
two subsystems. 
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Figure 2. Mn-projected d-DOS in Gai4Mn2NiQ (a) with Mn-Mn distance equal to 
3.19A and 5.19A and (b) 6.0A and 6.3A. The upper and lower panels represent the 
majority and minority spin DOS respectively. 

In order to understand the variation of exchange interaction among the TM 
impurity with the distance between the TM atoms, the concentration of Mn atoms 
was increased to 12.5%, equivalent to introducing 2 Mn atoms in the supercell. Self 
consistent calculations were carried out for two different magnetic configurations of Mn 
electrons in which the electrons are parallel or antiparallel corresponding to FM or AFM 
configuration. For all the geometries of the dopants, systems, as described in section 
II, with FM configuration of the Mn atoms were found to have lower energy. Since the 
Mn-d levels are responsible for observed half metallic behavior, a comparison of Mn-d 
DOS in various geometries is shown in figure El (a) for separations 3.19A and 5.2A and 



(b) for 6.0A and 6.3A. Here the half metalhc state is evident in all the cases. The 
TM-TM distance of 6.3A corresponds to single TM doping (6.25%). 

The Mn rf— DOS is broad for substitution at nn distance. In all the other cases 
the band is split indicating that the majority spin d-bands of the two Mn atoms at 
nn overlap. On increasing Mn-Mn distance, d-band splitting takes place implying a 
reduction in the interaction between the TM atoms. It may be noted that the minority 
spin conduction band overlaps with the majority spin band for nn substitution. A gap 
of ~ 0.5eV is present between the minority spin conduction band and majority spin 
band for Mn-Mn distance greater than nn. The minority spin valence band as well 
as conduction band is far apart from the ep thus retaining the highly spin polarized 
state also seen in the single Mn doped Gai^MnNiQ system. The down spin gap is 2eV 
for nn configuration and increases to 2.5eV at larger separations. On increasing the 
distance between the Mn atoms, splitting of d-level increases. This is consistent with 
the observation that in single impurity doping, TM-TM atom distance is 6.0A and 
splitting of the Mn-ci band is larger as seen in figure Efb). The magnetic moment at 
Mn-site does not depend on the distance between the dopant atoms and has a value 
3.34/iB for all the geometries as indicated in table [H 

Table 1. Magnetic moment (jUs) in various Mn-doped systems. Total„iom indicates 
mag. mom. /unit cell, dopantmom is the mag. mom. at dopant site and Nmom is the 
average mag. mom. at nn N sites. 



System 




Dopantmom 


^mom 


Gai^MnNiG 


4.00 


3.33 


0.001 


GaiiMn2NiQ 


8.00 


3.34 


0.001 


Mn-Mn = 3.19 A 








GaiiMn2NiQ 


8.00 


3.34 


0.005 


Mn-Mn = 4.5 A 








GaiiAIn2NiG 


8.00 


3.33 


0.006 


Mn-Mn = 5.19 A 








GaiiMn2Nie 


8.00 


3.34 


0.005 


Mn-Mn = 6.0 A 









The presence of localized moment influences the near neighbor N atoms within the 
GaN system, such that DOS of the nn N-atoms around the impurity atom becomes as 
shown in figure a). Due to the p — d interaction, induced states are seen on N atoms. 
The magnitude of the induced states is maximum at nn N atoms and decreases as the 
distance from the TM atom increases. This again indicates a localized nature of the 
TM states. Figure IHl^b) shows the spin charge density (SCD) in a plane containing three 
nn N atoms for single Mn doping. The plot shows that the SCD on the N atoms lying 
above the Mn atom is negative whereas it is positive for the N atom lying below the Mn 
atom. Average magnetic moment on the nn N atoms is positive as shown in table ^ for 
all the different geometries. 




Energy [eV] 



(a) 




(b) 



I <-0.0162 

[] -o.ooeo 

[] +0.0043 

n +0.0145 

n +0.024G 

[] >+0.0350 



Figure 3. (a) Variation in partial-DOS at 3 N-sites in Gai^MnNiQ. The solid lines: 
nn-N along c— axis. Dotted lincs:nn-N in plane. Dashed-dotted lines: next nn-N (b) 
SCD plot in a plane of three nn N atoms. 



3.2. Cr doped systems 

Electronic structure calculation for substitutional doping of Cr in the GaN system 
was also done and is analyzed in a similar fashion. For each Cr doped in the 32 
atom supercell, equivalent to 6.25% doping, there are five spin up d-states which are 
introduced in the GaN band gap. Since the Cr atom has 4s^3c?^ valence electrons, only 
three of the electron states out of the five ?> — d levels are occupied, thus creating two 
hole states per Cr substitution. Single Cr doping into the 32 atom supercell at cation 
site results in the Ci-d levels appearing in the band gap of the semiconductor host as 
seen in figure lU The Cr dtj levels split (figure EJ^a)), out of which two energy levels lie 
below the Fermi level {ep) and are occupied. The third level which is ~ 0.5eV above is 
unoccupied. One of the de level is occupied and the other lies just above the as seen 
in figure m^b). However, the hybridization of dt^ and d^ majority spin states is negligible 
and these levels are well separated as opposed to Mn doped case. As in the case of Mn 
doping, the Cr minority spin d-states are above the ep and so the impurity states at the 
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Figure 4. Cr-projected majority spin (a) and (b) bands in Gai^CrNiQ. 



ep are 100% spin polarized. 

In the Gai4Cr2NiQ system, figure El shows that even though the Cr atoms are 
substituted at nn sites there is a gap seen between the spht Ci-d levels, unlike in the 
nn Mn-doped system. The minority spin conduction band overlaps with the majority 
spin band for nn Cr case, as in nn Mn case. For all the geometries of the dopant atoms 
the system is half metallic. The band gap for the minority spin, in case of two Cr 
substitution at 3.19A is ~ 3eV which is larger than the corresponding Mn case which 
has a gap of ~ 2.0eV. When the TM-TM distance is 5.19A, in the Cr case the gap is 
~ 3.0eV whereas for Mn substitution it is ~ 2.5eV. The lowest energy configuration for 
Cr substitution occurs for Cr-Cr distance of 3.19A and for all the geometries studied 
the FM configuration of TM atoms has lower energy compared to AFM configuration. 
The magnetic moments at the Cr site in various geometries is as shown in table |21 and 
it is seen that the magnetic moment does not show much variation depending on the 
distance, indicating that the direct interaction between the Cr atoms is minimal. 



Table 2. Magnetic moment(/iB) in Cr-doped systems. Totalmom indicates mag. 
mom. /unit cell, dopant„iom is the mag. mom. at dopant site and f^mom is the average 
mag. mom. at nn N sites. 



System 




D op ant mom 




GaisCrNie 


3.00 


2.47 


-0.025 




6.00 


2.48 


-0.031 


Cr-Cr = 3.19 A 








Gai5Cr2Nie 


6.00 


2.47 


-0.022 


Cr-Cr = 4.5 A 








Gai5Cr2Nie 


6.00 


2.48 


-0.025 


Cr-Cr = 5.19 A 








Gai5Gr2Nie 


6.00 


2.48 


-0.026 


Cr-Cr = 6.0 A 
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Figure 5. Cr-projected d-DOS in Gai4Cr2-/Vi6 for Cr-Cr separation of (a) 3.19A and 
5.19A and (b) 6.0A and 6.3A. 




Figure 6. SCD plot in a plane of three nn N atoms around the TM atom in 
GaiiCrNiQ. 



Magnetic moments observed at the various dopant sites from our calculation are 
shown in the table El The total magnetic moment per unit cell per Mn atom is A^b and 



the average magnetic moment on the nn N atoms in case of Mn-doping is parallel to the 
Mn-moment. This can be understood as penetration of the spin-polarized Mn states 
to the neighboring host which does not have any of its own states in the gap region. 
The magnetic moment per unit cell per Cr atom is S/zs. The average magnetic moment 
on the nn N atom is anti-parallel to Cr- moment. The difference in the orientation of 
the average magnetic moment on nn N atoms of Mn and Cr is due to the difference 
in the p DOS of the nn N along the z axis (figure not shown here) compared to the 
nn N atoms lying in the xy plane above the TM atoms. There is not much variation 
of magnetic moment with increase in Cr-Cr distance and shows a similar trend as Mn 
doped systems. The magnitude of average nn-N magnetic moment is greater in case 
of Cr substitution, which contribute one less electron to the hybridized valcen band. 
SCD in figure El on all of the nn-N atoms of the single Cr (only 3 nn N atoms shown in 
figure El) doped system is negative thus showing that the TM atom and the nn-N are 
anti-ferromagnetically coupled. 

3.3. Estimation ofTc 

We have predicted the for the DMS based on Gai6Mn2Ni6 and Gai6Cr2Ni6 considering 
the mean field approximation. Figure shows the AE = Eafm — Epu for the Mn/Cr 
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Figure 7. /S.E — Eafm — Efm for Mn and Cr doping. Inset shows the mean field 
Tc variation with distance between dopants. 

doped systems, where Eafm is the total energy for the antiferromagnetic (AFM) 
configuration and Efm is the total energy for the ferromagnetic (FM) configuration. 
Observed variation AE vs distance for Mn substituted DMS agrees with the one reported 
by Sanyal JHl ■ AE is a measure of the exchange interaction in the system. Highest AE 
is seen for the case where the TM atoms are substituted as near neighbors, signifying 
larger overlap of the magnetic impurity orbitals. For Mn doping at nn AE ~ 0.14eV, 
this compares well with the value calculated for dimer substitution by Uspenskii et. 
al ^21 which compared qualitatively with the high = 940K measured P|. As for 
the identical Cr case AE ~ Q.lQeV also compares well with the observed but is 




lower than the 900K observed by Liu and co-workers j^. From figure [7| it also emerges 
that the exchange interaction decreases sharply as the distance between the TM atoms 
increases. Thus the exchange interaction is short range and could be interpreted as the 
double exchange mechanism. 

4. Summary and conclusions 

We have analyzed the electronic structure of GaN doped with TM Mn and Cr with 
6.25% and 12.5% doping for various possible geometries to replicate the situation where 
the TM atoms would appear either to cluster or be separated. The self consistent 
FP-LAPW calculations predict half metallic state for 6.25% as well as 12.5% doping. 
Comparing the total energies of the FM and AFM configurations for 12.5% doping, 
the FM state is found to be lower in energy and is predicted to be the preferred state. 
On-site magnetic moment at the TM site shows insignificant variation with distance 
between the dopants. The near neighbor N atoms contribute to the states in energy gap 
of the semiconductor due to the infiuence of the TM atoms. Average magnetic moment 
at nn N site is parallel to the Mn magnetic moment where as it is anti-parallel to the Cr 
atoms. We observe that both the systems with nn substitution of Mn/Cr atom would 
show high Tc. The energy gap between the minority spin band in Mn is ~ leV lower 
than in Cr doped system and we think this could be an important factor in determining 
a suitable system. But since the magnetic moment at Mn site is higher than Cr, it 
would be of interest to study mixed systems of Mn and Cr to incorporate the salient 
features of both TM atoms. 
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